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Abstract 

The anyonic quantum walk is a dynamical model describing a single anyon propagat- 
ing along a chain of stationary anyons and interacting via mutual braiding statistics. We 
review the recent results on the effects of braiding statistics in anyonic quantum walks 
in quasi-one dimensional ladder geometries. For anyons which correspond to spin-1/2 
irreps of the quantum groups SU{2)k, the non-Abelian species (1 < fc < oo) gives rise 
to entanglement between the walker and topological degrees of freedom which is quan- 
tified by quantum link invariants over the trajectories of the walk. The decoherence is 
strong enough to reduce the walk on the infinite ladder to classical like behaviour. We 
also present numerical results on mixing times of SU {2)2 or Ising model anyon walks 
on cyclic graphs. Finally, the possible experimental simulation of the anyonic quantum 
walk in Fractional Quantum Hall systems is discussed. 

1 Introduction 

In two spatial dimensions, the exchange statistics of identical particles must be generalized 
to braiding statistics pLj. The particles which obey braiding statistics are called anyons, and 
they are thought to exist as quasiparticle excitations in Fractional Quantum Hall systems [2] 
or as edge modes of semiconductor nanowires [3] . Much of the interest in anyons is motivated 
by the realization that particles with non-Abelian braiding statistics can be used to perform 
quantum computation in a decoherence-free manner [1]. In such schemes the information is 
encoded in the collective many-body state of a system in a topological phase of matter. 

Information stored in the collective state of anyons can be manipulated by either fusing or 
braiding anyons. Recently there has been some research activity on the equilibrium properties 
of systems of anyons interacting via fusion [SJ El E]- Here we discuss the non-equilibrium 
properties of anyons which interact via braiding. The dynamics of anyons can be simulated by 
a discrete-time quantum walk, where the motion of a moving anyon is coupled to an external 
degree of freedom (DOF) called the quantum coin. The moving anyon, called the walker, 
braids with an array of stationary anyons arranged on a line. We concentrate on the effects 
of braiding on the propagation of the moving anyon along the lattice. 
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Figure 1: The anyonic quantum walk model. Here the coin modes and 1 label the position 
of the walker (red cross) on either upper or lower edge, respectively. The locations of the 
contacts between the edges are labelled by spatial sites s. The coin flip matrix U corresponds 
to a scattering matrix at the edge contacts. The action of braiding between the walker anyon 
and the stationary anyon (grey cross) between sites s and s + 1 is represented by the generator 



If the anyons are non-Abelian, the mobile anyon and the stationary anyons possess a 
collective Hilbert space Hp, which grows exponentially in the number of anyons, and which 
allows for interactions between them. In fact this coupling induces decoherence on spatial 
DOF of the mobile anyon, and we study the effects of this coupling in both finite and infinite 
lattices. Differences between Abelian and non-Abelian anyons are pointed out, and possible 
implications for distinguishing between these two types of anyons experimentally are discussed. 

Phase factors in two-particle walks, corresponding to braiding phases of Abelian anyons 
were studied by Berry and Wang [12j and a similar scheme was implemented experimentally 
with integrated photonics by Sansoni et al [13]. Also there has been some related work on 
using quantum walks to simulate different phases of topological insulators [H [9l [TOl [11] . In 
such systems, a discrete walk models coarse grained single particle continuous dynamics of a 
system possessing bound states at the boundaries of different symmetry protected topological 
phases. These bound states are robust to small perturbations of the system which respect the 
symmetry (e.g. time reversal symmetry) and can be probed by varying the step dynamics in 
the quantum walk. In contrast our model involves multiple particles with topological symmetry 
which is robust to arbitrary local perturbations. 

2 Anyonic quantum walk 

The dynamical behaviour of anyons can be studied by decimating the time evolution into 
discrete steps of infinitesimal length. Here we consider a single particle (walker) hopping 
on a one-dimensional lattice of spatial sites. Between each discrete time step, the walker is 
allowed to move only to its neighbouring sites on the left and right. To study the effects of 
anyonic interactions in the dynamics, we place an array of stationary anyons on the dual lattice 
between the spatial sites as shown in Fig. [T] The mobile anyon and the stationary anyons 
are assumed to stay far enough apart from each other such that the interactions are mediated 
exclusively by braiding statistics. Although the spatial lattice extends in one dimension only, 
it is necessary to think of the system as two-dimensional, such that the mobile anyon can 
circumnavigate the stationary anyons without coming into contact with them. 

The braiding interactions of anyons are purely topological. The action on the wave function 
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does not depend on the exact path of the walker, only on the initial and final points For 
simplicity of analysis we consider only counterclockwise braids, meaning the system is chiral. 
Numerical results on continuous-time anyonic walks show that the essential features of the 
walk do not change for non-chiral walks [211 E] • 

The formalism of anyonic quantum walks [T5] is a generalization of the discrete-time 
Hadamard walk on an infinite one-dimensional lattice as defined by Ambainis et al [16]. In 
the anyonic quantum walk model the Hadamard walk is changed to include the effects of 
the braiding interactions. In addition to the spatial and coin degrees of freedom, one has 
to accommodate for the fusion Hilbert space of (non-Abelian) anyons, and the total Hilbert 
space is thus T-L = Tis ® Tip ® 'He- The space T-Ls of spatial sites is spanned by the basis 
vectors Tip is described below, and the two-dimensional coin space Tic is spanned by 

{|0), |l)}, where corresponds to going left and 1 to going right. The system is initially in the 
product state |^'(0)) = |so)g|$o)p|co)Q- In the initial configuration, the walker is localized 
at some initial site Sq and the coin state is chosen to be |co) = |0). The initial state of the 
fusion space |$o) is the vacuum state such that the anyons taking part in the walk are created 
in particle-antiparticle pairs. If the particle-antiparticle pairs reside on neighboring islands, 
the tracing scheme is determined to be the plat closure. If the antiparticles are thought to 
be ancillary particles which are dragged out of the system for the duration of the walk, the 
tracing scheme is Markov. 

The time evolution between infinitesimal time steps is given by the operator |^E'(t + 1)) = 
Vl^|\E'(t)) which consists of the coin fiip operator and a subsequent conditional braiding op- 
erator W=TU. The coin fiip operator acts like the Hadamard gate on the coin DOF: 

U = Is Ip ® , and the conditional braiding operator moves the walker to the 

adjacent lattice sites, conditioned on the coin state and simultaneously braiding with the 
stationary anyon between the lattice sites (see Fig. [T]): 

T{\s)\f)\c)) = \s- l>6,_i|/>|0>(0|c> + \s + l>6,|/>|l>(l|c>. (1) 

The wave function at time t is obtained by applying the single step operator repeatedly on 
the initial state = iy*|^(0)). 



2.1 Representations of the braid group 

The generators of the braid group bs E B^^ act in the fusion space T-Lp only. The space Tip is 
spanned by vectors {\gi, g2, . . .)} where the quantum numbers gj label possible intermediate 
fusion outcomes of identical anyons. The representations of bg depend on the anyon model, 
which is defined by the set of all charge types and the fusion rules of the charges. Here we 
give the representations of bg for two different particle types, the a charges of the so called 
Ising anyon model, and the spin-1/2 irreducible representations (irreps) of the quantum group 
SU (2)2. The fusion spaces of these particles have a very convenient tensor product structure in 
terms of qubits. Every even intermediate fusion charge of an A^^-anyon state in the standard 
basis is constrained to be g2j=o' (or g2j=\) for the Ising and SU{2)2 models respectively, 

""^ There could be dynamical and or geometric phases arising from e.g. Aharanov-Bohm interactions but we 
neglect them in the following. Provided the interactions that give rise to these phases are quenched in time, 
they do not qualitatively affect the transport behaviour of non-Abelian anyons. If they are random in space 
they can induce localization in Abelian anyons j21j . 



3 



and every odd charge is constrained to two values g2j+i^{l,ip} (or g2j+i&{0, 1}), forming a 
single qubit. Requiring that the total charge is vacuum, the fusion space then consists of 
m = Na/^ — 1 qubits, T-Lp ~ (C^)*^"^. The braid generators of Ising and SU {2)2 anyons thus 
admit a simple structure where the non-trivial matrices acting on the fusion space are at most 
4-by-4 [15]. The even and odd generators can be expressed as 

hi = R®J=2h, b2 = B(g)JL2h, h = P(g)p„ 

b2k = ®j=i^l2 ® B 12, &2fc+i = ®j=Jl2 ® P ®f=fc+2 h; l<k<m, (2) 

For the a anyons in the Ising model the non-trivial matrices are 

i? = P.,. = e-i (^J B = ^(^^ll^ ^ = e-^tdiag(l,z,z,l). (3) 

The R symbols of the SU {2)2 model, describing the braiding of two particles, are obtained 

from those of the Ising model by the substitution R\ ^ = —iR^^a where the star denotes 

2 ' 2 

complex conjugation. The F symbols, which describe changing the order of fusion for three 

particles, are identical for the two models, explicitly F = Fi^ii = ^ ^j- Thus up 

to phase factors and complex conjugation, the Ising and SU{2)2 models act similarly under 
braiding of spin-1/2 irreps with charge correspondence {1, cr, ■?/'} = {0, ^,1}. Consequently the 
braid generators for the SU{2)2 anyons are given by Eq. ^ and substituting the matrices in 
Eq. (g with R^ iR*, B iB* and P ^ iP*. In most cases the difference between these 
two models is not distinguishable. 



2.2 Probability distribution and quantum link invariants 

The main interest in the anyonic quantum walk is in how the occupation probability of the 
mobile anyon on the lattice sites evolves, if it is initially localized at some initial site sq. The 
probability to be at site s at time step t is given by the diagonal values of the reduced density 
matrix: 

pis,t) = (psW),,, = {s\TTFTTc{W'\^iO)){^iO)\iW^Y)\s). (4) 

In a t-step walk, the total transition amplitude between sites sq and s can be given as a sum 
over all the possible ways to end up on site s after taking t steps: 

P(^^^) = I E (-1)^^'^+^^''^ Tryid,d',%) (5) 

where paths are labelled by vectors a = {ai, a2, . . . , at} with aj=0 {aj=l) if the walker went 
left (right) at time step j. The quantity z{a) = X]j=i '^i'^i+i counts the number of subsequent 
right-moves in the path a, and the anyonic term is given by 3^(a, a',$o) = 'SQ:|$o)($o|'Si/ 
with the braid word Bs given as a product of generators corresponding to path a: Bs = 
nt=o ^so+2{x;*"i ai)-(t-r)' "^^^ notation a.a' s means that the summation is done only over 

paths which satisfy X]j=i — Sj=i ^'j ~ ''"2'"'"* ^^"^ ~ constraints which arise from 
taking the diagonal elements of the spatial DOF and tracing over the coin. The trace over 
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the anyonic DOFs for spin-1/2 irreps of SU{2)k can be expressed in terms of link invariant 
polynomials, as shown by Aharonov et al [17]: 

('_o3/4\«,(L(a,a')) 

Tr3^(a,a',<|.o) = (<fo|4S^|$o> = ^ \L_, V^l{.,^')(?) (6) 



(7) 



where Vi^q) is the Jones polynomial of a link L with parameter q = e**+2, {|L(a, a'))(A) is 

1 //I — ^TT 

the Kauffman bracket with parameter A = q~ ' = e'^^^+^\ and the quantum dimension is 
d = —A? — A~'^ = 2 cos By choosing the initial state $o suitably, namely as a product 
of vacuum pairs with one member of each pair taking part in braiding, the link L{a,a') is the 
Markov closure of the braid word B-,Ba- 



2.3 Infinite chains 

When there are no boundaries in the system, the walker can extend to infinite lengths on 
the left and right. If the walker is initially confined in some finite region, locality of the shift 
operators implies that there is a maximum distance that the walker can propagate in a finite 
time t. The main question is that if the walker is initially localized at some initial site sq, 
what is the expected distance of the walker from so at time step t? The square of the expected 
distance is quantified by the variance 

a\t) = {{s- sof) = t)s^ - 2so t)s + si (8) 

s s 

One can identify different qualitative behaviour of walks by the functional dependence of the 
variance on the number of time steps t. If cr^(t) = K2t'^ + Kit + Ko the mean velocity of the 
walker is linearly proportional to the maximum allowed velocity, and the behaviour is called 
ballistic. This is the behaviour of the Hadamard walk without decoherence. Random walks 
and quantum walks with decoherence scale as cr^(t) = K^t + Kq and propagation is called 
diffusive. Localization implies that cr^(t) < Kq, but the exact definition of localization can 
vary [22] . 

For Abelian anyons, the walker acquires a global phase e'^'^^'^ from each braid with the 
stationary anyons. The total phases for bra- and ket-evolution of the density matrix are e**'''/^ 
and e"**'''/^ respectively, so that the total effect is trivial and the braiding statistics plays no 
role in the time evolution. Non-trivial effects do arise for Abelian anyons if one introduces 
clockwise braiding and extends the coin space to four dimensions. In this case the behaviour is 
ballistic for all values of the exchange angle (f, but the coefficient of the leading term depends 
on [I5] . 

The anyonic quantum walk of SU {2)2 anyons was analysed in Ref. [18|. For the value of 
the parameter q = i, the Jones polynomial can be evaluated efficiently and the probability 
distribution can be written as 

r ^_iy{a)+zia')+r{a,a') L proper , , 

^ ' ' J-^ 1 L not proper ^ ' 

where properness and T{d,a') are properties of the link corresponding to the path (a, a'). 
Numerical results show that the walk is diffusive, and the variance was shown to be linearly 
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dependent on t at the asymptotic limit t — )■ oo. These results hold equivalently for SU {2)2 
and Ising anyons. 

The classical diffusion like behaviour of the non-Abelian Ising anyons is not special to 
that model but in fact is generic for a large class of non-Abelian anyons. To show this we 
investigated quantum walk behaviour of spin— 1/2 irreps of all SU{2)k models. For k > 2 and 
finite these models have non-Abelian anyons. An obstacle to the analysis, however, is that for 
higher values of fc, the evaluation of the Jones polynomial is exponentially hard in the number 
of time steps, and the algorithm that was used for k=2 does not work. Numerical results for 
i:=10 show that the distribution looks classical for small values of k and approaches quantum 
distribution when k increases [15]. Results for higher number of time steps can be obtained by 
adopting the so called W'^ model [19] where the walker evolves according to the time evolution 
operator W for two time steps, after which the fusion and coin DOF are traced out. A single 
step of the walk is described by the completely positive trace preserving map 

ps(t + 1) = S{ps{t)) = J2 Efc Ps{t) El (10) 

f,c 

where Efc are Kraus generators, labelled by the basis states of the fusion and coin DOF. 
Numerical results show that the walk is also diffusive for all values of k. Approximative 
expression for SU (2) 2 /Ising anyons for large t can be obtained by approximating the Kraus 
generators as circulant matrices. The result is 

1 / t 

which is the probability distribution of the classical random walk. 

The behaviour of non-Abelian walks can be explained by decoherence effects in quantum 
systems. The spatial and coin DOFs become entangled with fusion DOF, and the state is 
not pure anymore in the space Hs ® ^c- This can be seen by considering the representations 
of the generators of the braid group as given in Eq. ([2]). When the walker propagates along 
the lattice, it braids with a growing number of stationary anyons, and the fusion qubits are 
multiplied by matrices R, B and P. Consider the state of two previously untangled fusion 

qubits i and j: l^i) = (alo) + /3|l)).(a|0) + ., with a = {Fhi)o and /3 = {Fhi)i- 

' ' VI/ 222 222 

These could represent one encoded qubit at the current position of the walker on the fusion 
tree and the other among adjacent anyons not yet braided around by the walker, represented 
in a basis where the two fusion qubits belong to separated fusion trees with three charges 
each and total charge for all six equal to zero [20]. In the sector of | total charge for each 
tree, the single intermediate charge of each tree defines a qubit. The walker becomes linked 
nontrivially with a stationary anyon if it braids twice with the stationary anyon. This operation 
is represented by the square of their mutual braid generator, \ip2} = The amount of 

entanglement between the fusion qubits generated by the braid 6^ can be quantified by the 
linear entropy of the reduced density matrix of the qubit i: Q{pi) = 2(1 — Tr(p^)/Tr(pi)^) 
where pi = TTj[Pij\4)2){'4'2\]- Here Pij is a projector onto the two qubit subspace, necessary 
since the braiding can couple to a fifth leakage state with total charge in each tree equal to |, 
and the linear entropy normalises by the trace of the reduced state pi. We find that the linear 
entropy is maximal Q{pi) = 0.938 for the level fc = 3 of Chern-Simons theory, and decreases 
as a function of fc, with an asymptotic behaviour Q{pi) = — 0(l/fc^). This is consistent 
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with the fact that the representations of SU{2)k approach those of SU (2) at the hmit of large 
k, and the the action of two-fold exchange is trivial in SU{2) (no entanglement generated). 
Furthermore, the linear entropy is zero for k = 1, which corresponds to Abelian anyons by 
fusion rules. The fusion space dimension of Abelian anyons is one and no entanglement is 
therefore generated by braiding. 

All the results discussed here hold for uniform configurations of anyons, ie. there is exactly 
one anyon with multiplicity 1 between each spatial site. More generally, one could define that 
each island between the lattice sites can host more than one anyon. The transport properties 
of anyons in such random charge configurations were considered in Ref. [21] with dynamical lo- 
calization phenomena in mind. This might be relevant for experimental simulations of anyonic 
quantum walks since the quasiparticle occupations might be hard to control in experimental 
conditions. 

3 Finite chains 

Finite lattices can be modelled by truncating the graph on the left and right and introducing 
boundaries at these locations. This requires a choice of boundary conditions for the time 
evolution operator W. Possible choices for unitary walks are periodic and reflective boundary 
conditions, and for non-unitary walks, absorbing boundary conditions. More generally, one 
could think of the boundary as an interface between the quantum system and its environment, 
such that the walker interacts with the environment at the boundary, and this interaction 
changes the state of the walker in some way. 

The reflective and absorbing boundary conditions can be defined in a straightforward way 
by introducing ancillary states at the boundaries [22] • The stationary anyons reside only in 
the bulk, so the braiding is trivial when the walker moves between the bulk and the boundary 
states. For periodic boundary conditions the situation is slightly different, since the transitions 
between lattice sites 1 and N must be accompanied by changes of labels of the fusion charge. 
For example, if the walker moves from site to site 1, the fusion label changes to 1 
and all other labels are added by one. This operation is represented by the braid word 
Bna^i = n^i""'^ bk and the inverse operation is given by Bi^Na = -SjvA^i- 

The anyonic quantum walk on a finite chain can be simulated numerically by constructing 
the initial state |\E'(0)) = |so)|$o)|co) and applying the time evolution operator t times, 
|\E'(t)) = H^*|\&(0)). The probability distribution at time step t is given by the diagonal values 
of the reduced density matrix of the position space, see Eq. Q, and the braid generators bs for 
Ising and SU{2)k anyons can be constructed using Eq. It is known that on finite graphs, 
the quantum walk does not converge into any asymptotic distribution (unless the system is 
initially in an eigenstate of the time evolution operator, in which case the walk is trivial). 
However, the time-averaged probability distribution, 

I * 

Pis,t) = —J2pis,T) (12) 

r=0 

has been shown to converge [23] , but the asymptotic distribution is not necessarily the uniform 
distribution. As the classical random walk is a Markov process, any information on the 
initial state is lost after a certain number of time steps, and the probability distribution 
becomes uniform on the nodes in the infinite time limit. One can therefore compare the 
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convergence of the random walk and the time-averaged quantum walk to their respective 
asymptotic distributions. 

3.1 Mixing time 

The convergence to the asymptotic distribution can be quantified by different measures. One 
such measure is the e-mixing time defined as 

M, = min{r| IIA < e, Vt > T} (13) 

where Dt is the probability distribution at time step t, the asymptotic distribution is vr = 
lim Dti and the distance measure between two probability distributions is the total variation 

distance defined as \\Di — -D2II = Xli ~ D2{i)\- The mixing time of random walks is 

shown to converge as ~ 0(A^^ log 1/e). Aharonov et al. showed that the quantum walk 
converges as ~ 0{N \ogN 1/ e^), which provides a quadratic speedup in number of sites N , 
but an exponential slowdown in the parameter e. The convergence with respect to e can be 
improved by amplification, and Richter has shown that a suitably chosen amplification scheme 
can enhance the convergence to 0(A^logl/e) |24 ]. 

As the asymptotic distribution is not currently known for anyonic walks, an operational 
measure of convergence is defined as the total variation distance to the final distribution at 
time step T: 

D(t,r) = ||^(M)-pM1|| (14) 

such that lower values of D mean faster convergence. The total variation distance for Ising 
anyons is plotted in Fig. |2]with the corresponding plots of the random walk and the Hadamard 
walk (also the average distribution of the random walk is plotted for comparison). For a small 
number of sites, the Ising walk converges in a similar fashion as the Hadamard walk. The total 
variation distance is initially similar to the classical random walk (RW), and fiuctuates between 
RW and RWavg at later time steps (mixing time for RW is calculated using the instantaneous 
probability distribution and RWavg is calculated using the time-averaged distribution, as with 
quantum walks). For a larger number of sites, the convergence of the Ising walk is still similar 
to the Hadamard walk. Both walks converge faster than the RW initially, but in the course of 
time the total variation distance of both the Ising and Hadamard walks settles between RW 
and RWavg- For a large number of sites, the convergence of the Ising walk is quite smooth, 
whereas the convergence of the Hadamard walk is not uniform. 

It is curious to note that the mixing properties of the non-Abelian anyonic walk are quite 
similar to those of the quantum walk. This is in contrast to the results on infinite chains as 
discussed in the previous section, where the propagation of the non-Abelian anyon shows sim- 
ilar properties to the classical random walk. Analysis of the probability distributions in finite 
chains shows that both the anyonic and non-anyonic quantum walks evolve in a disordered 
manner, and the probability distributions between individual time steps can be very differ- 
ent, whereas the probability distributions of the classical random walk are similar. The Ising 
and Hadamard walks sample from a wider range of probability distributions, and therefore 
the average distribution converges slower than RW. However, the Ising and Hadamard walks 
converge faster than RWavg- 
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Figure 2: Total variation distance D{t,T) between probability distributions at time steps t 
and T. a) 5 sites, b) 21 sites. Black represents the classical random walk and green represents 
the time-averaged classical random walk. Blue is the time-averaged Hadamard walk and red 
is the time-averaged Ising walk. Total number of time steps was 1000, periodic boundary 
conditions were used, and the initial fusion state |$o) was chosen such that tracing scheme is 
plat closure. 

4 Proposal for experimental setup 

A possible experimental platform for simulating anyonic walks could be in Fractional Quan- 
tum Hall samples. The quasiparticles at the z/ = 5/2 plateau have been proposed to have 
non-Abelian anyonic statistics corresponding to Ising anyons. The experiments probing quasi- 
particle statistics can be done using a Fabry-Perot type interferometer on top of a FQHE layer 
[251 [26] . An interference pattern is induced on the longitudinal conductivity (T^x as the side 
gate voltage Vs is varied. 

The statistical phases of anyons are generally very hard to detect. The interference pat- 
tern of a^x is dominated by Aharonov-Bohm oscillations, the period of which depends on the 
effective electric charge e* of the quasiparticles in the medium. The effective charge of quasi- 
particles at the z/ = 5/2 plateau is e* = e/4 in all proposed models, regardless of whether 
the statistics is Abelian or non-Abelian, so the period of Aharonov-Bohm oscillations does 
not distinguish between the statistics. It was proposed however that if the quasiparticles are 
non-Abelian and if there is an odd number of them between the tunneling point contacts of 
the interferometer, the Aharonov-Bohm oscillations disappear [271 EB] • The disappearance of 
the oscillations was confirmed experimentally by Willett et al [29], but instead of complete 
suppression of the oscillations there seem to be persistent e/2 oscillations when the number of 
quasiparticles is odd. The e/2 oscillations could result from tunneling of Abelian quasiparti- 
cles with effective charge e* = e/2, or non-Abelian quasiparticles tunnelling more than once. 
Analysis by Bishara et al. shows that the first option is more plausible [5U] . 

4.1 Setup for measurements 

We propose that an experimental setup which implements the anyonic quantum walk scheme 
could also be used to probe the statistical properties of quasiparticles. Such a scheme relies 
on transport measurements of edge mode quasiparticles in a multipoint contact version of 
the Fabry-Perot type interferometer [18j. Instead of measuring the interference pattern of 
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conductivity when the side gate voltage is varied, the conductivity could be measured in 
equilibrium while keeping the side gate voltage constant. We suggest that the quasiparticle 
statistics has an experimental signature on the conductivity along the point contact array, and 
that this signature becomes clearer when the number of tunneling point contacts increases. 

To measure the probability distribution of the quantum walk directly would require in- 
troducing contacts at each location corresponding to spatial sites of the quantum walk. In a 
simplified scheme, very similar to that used in two-point contact interferometry, the conduc- 
tivity would only measured at the ends of the multipoint contact ladder. Quasip articles are 
ejected to the array from the middle or from the other end via an external contact. Changing 
the voltage by the chemical potential of the quasiparticles creates pairs of quasiparticles on 
the edges. Thus, their wavepackets are initially localized at the location of the contact, and 
they start to perform a quantum walk across the ladder. To measure the conductivity between 
two points, one has to introduce another contact at some other point in the ladder. If the 
measurement destroys the quasiparticle, this corresponds to introducing an absorbing bound- 
ary to the quantum walk. The conductivity between the ejection point and the measurement 
point is proportional to the probability to reach the absorbing boundary when the walker 
starts localized at the ejection point. 

4.2 Exit probability 

To illustrate how statistics can affect the observed conductivity, consider a scheme where 
quasiparticles are ejected from the middle of the ladder, and the sum of conductivities at the 
ends of the ladder is measured. In principle the conductivity is related to the probability to 
propagate from the initial site to the edges of the array. Starting from the initial site sq at 
t=0, this probability is called the accumulated exit probability defined by 

t 

Pex(t) = 5^p(0,r)+p(iV + l,r), (15) 

T=0 

where the sites s = {0, + 1} are absorbing boundary sites. The exit probability is plotted 
in Fig. |3] as a function of time for the classical random walk, the Hadamard walk and the 
Ising anyon walk. At short time scales the Ising walk and RW behave quite similarly. The 
probability to observe the particle is higher for QW, as is expected for ballistic propagation 
of the walker. At long time scales the probability to observe the RW particle is the highest 
and it approaches 1. The probability to observe the Ising anyon is the lowest. Both the Ising 
walk and the QW probabilities seem to approach a constant < 1, which is in agreement with 
earlier results which show that the asymptotic exit probability does not approach unity [16j . 
It should also be noted that these effects become clearer when the number of sites A^ grows. 

A few problems with these proposals can be immediately identified. The results presented 
here use the Hadamard coin which corresponds to strong backscattering limit at the quantum 
point contacts, which is known to be problematic in observation of statistics. However, the 
main results shouldn't be too sensitive to the choice of coin because even for small coin 
rotation angles, there will entanglement generated for non-Abelian anyons which will affect 
the transport behavior in a similar manner though over a longer time scale. Also, it might 
be hard to protect the mobile quasiparticles from spatial decoherence, such that they would 
propagate coherently across several point contacts across the sample. If this were the case 
then Abelian and non-Abelian anyons would give similar signatures. Finally, the quasiparticle 
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Time step t Time step t 



Figure 3: Accumulated exit probability as a function of time. The number of sites is a) = 12 
and b) = 28. 

occupations on the islands and gate voltages might fluctuate between experimental runs. In 
this case the model would have to take the effects of disorder into account [2T]. 

5 Discussion and outlook 

The effects of particle statistics in the one-dimensional anyonic quantum walk were analysed, 
and the results on the infinite line indicate that the walker propagates diffusively, showing 
similarities with classical random walks and quantum walks with decoherence. In finite chains, 
the effects are less pronounced. The mixing properties of the average probability distribution 
are close to those of the Hadamard quantum walk, and the probability to be absorbed at 
the boundaries is between those of RW and QW. Experimental realizations using an array of 
quantum point contacts was also discussed. 

One interesting generalization of the models presented here is obtained by considering 
non-uniform charge configurations [21]. Some work is also being done with continuous-time 
anyonic quantum walks using matrix product states and Hubbard-type Hamiltonian models. 
Other natural generalizations include multiparticle walks with more than one mobile anyon 
and anyonic quantum walks in two-dimensional lattices. 
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